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We investigate the dynamics of polymer translocation through a nanopore using two-dimensional 
Langevin dynamics simulations. In the absence of external driving force, we consider a polymer 
which is initially placed in the middle of the pore and study the escape time r e required for the 
polymer to completely exit the pore on either side. The distribution of the escape times is wide and 
has a long tail. We find that r e scales with the chain length N as r e ~ N 1+2 ", where v is the Flory 
exponent. For driven translocation, we concentrate on the influence of the friction coefficient £, the 
driving force E and the length of the chain N on the translocation time r, which is defind as the 
time duration between the first monomer entering the pore and the last monomer leaving the pore. 
For strong driving forces, the distribution of translocation times is symmetric and narrow without 
a long tail and r ~ E~ x . The influence of £ depends on the ratio between the driving and frictional 
forces. For intermediate £, we find a crossover scaling for r with N from r ~ N 2u for relatively 
short chains to r ~ N 1+v for longer chains. However, for higher £, only r ~ N is observed even 
for short chains, and there is no crossover behavior. This result can be explained by the fact that 
increasing £ increases the Rouse relaxation time of the chain, in which case even relatively short 
chains have no time to relax during translocation. Our results are in good agreement with previous 
simulations based on the fluctuating bond lattice model of polymers at intermediate friction values, 
but reveals additional features of dependency on friction. 



I. INTRODUCTION 



The transport of a polymer through a nanopore is as- 
sociated with an energy barrier arising from loss of con- 
figurational entropy due to the geometric constriction. 
Such processes are commonly observed in biology, such 
as DNA and RNA translocation across nuclear pores, 
protein transport through membrane channels, and virus 
injection H 0. Moreover, translocation processes 
have various potential technological applications, such as 
rapid DNA sequencing HJali gene therapy and con- 
trolled drug delivery, etc In addition, the translo- 
cation dynamics is also a challenging topic in polymer 
physics. Particularly, the scaling of translocation time 
r with the chain length N is an important measure of 
the underlying dynamics of polymer translocation. As 
a result, recently a considerable number of experimen- 
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have been carried out. 



In order to overcome a large entropic barrier in poly- 
mer translocation and to speed up the translocation, an 
external driving force is needed in experiments, such as 
an electric field, chemical potential difference, or selective 
adsorption on one side of the membrane. For example, 
in 1996, Kasianowicz et al. reported an elegant ex- 
periment where an electric field can drive single-stranded 
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DNA and RNA molecules through the a-hemolysin chan- 
nel of inside diameter 2 nm and that the passage of each 
molecule is signaled by the blockade in the channel cur- 
rent. In fact, the translocation process includes two es- 
sential steps. First, one end of the polymer enters the 
pore directed by diffusion and by the action of an elec- 
tric field near the pore. The experimental results show 
that the ability of the polymer to enter the nanopore 
depends linearly on polymer concentration and exponen- 
tially on the applied voltage . Second, the poly- 
mer is translocated from one side of the membrane to the 
other, driven by the electric field. In the experiment for 
a-hemolysin, the linear behavior of r with the N is ob- 
served I n addition, an inverse linear dependence 
and an inverse quadratic dependence of the translocation 
time on applied voltage are observed for different exper- 
iments p7L fTT| . 

To overcome the limited voltage range that can be ap- 
plied across a biological pore and the difficulty in analyz- 
ing the current variations because the shot noise is com- 
parable to the expected signal, recently Li et al. fl4l Il5j 
showed that a solid-state nanopore could also be used 
for similar experiments. Most Recently, Storm et al. 01 
carried out a set of experiments on double-stranded DNA 
molecules with various lengths that translocate through 
a solid-state nanopore 16]. Surprisingly, a power-law 
scaling of the most probable translocation time with the 
polymer length, with an exponent 1.27 was observed, in 
contrast to the linear behavior observed for all experi- 
ments on a-hemolysin channel. 

Inspired by the experiments F?l Hll Pill a number of 
recent theories JlMlJl|j2ClJ^ IH H 
H3, IH, Ha Ha H3, IH, HE HI, HE Hl| have been devel- 
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oped for the dynamics of polymer translocation. Even 
for the field free case, polymer translocation remains a 
challenging problem. To this end, Sung and Park [Tfll | 
and Muthukumar p2^ | considered equilibrium entropy of 
the polymer as a function of the position of the polymer 
through the nanopore. Standard Kramer analysis of dif- 
fusion through this entropic barrier yields a scaling pre- 
diction of the translocation time r ~ N 2 for long chains. 
However, as Chuang et al. [13 noted, this quadratic scal- 
ing behavior cannot be correct for a self-avoiding poly- 
mer. The reason is that the translocation time is shorter 
than the equilibration time of a self-avoiding polymer, 
Tequii ~ N 1+2u , where v is the Flory exponent [461 l47j . 
thus rendering the concept of equilibrium entropy and 
the ensuing entropic barrier inappropriate for the study 
of translocation dynamics. Chuang et al. 32] performed 
numerical simulations with Rouse dynamics for a two- 
dimensional (2D) lattice model to study the transloca- 
tion for both phantom and self- avoiding polymers. They 
decoupled the translocation dynamics from the diffusion 
dynamics outside the pore by imposing the artificial re- 
striction that the first monomer, which is initially placed 
in the pore, is never allowed to cross back out of the pore. 
Their results show that for large N, translocation time 
scales approximately in the same manner as equilibration 
time. 

For forced translocation, Sung and Park 0] and 
Muthukumar [22j suggested a linear dependence r ~ 
TV under a strong field. This is in agreement with 
some experimental results 0, 1 1 lj for polymer translo- 
cation through a-hemolysin channel. In addition, dy- 
namic Monte Carlo (MC) simulation using Gaussian 
chain model [37j and Langevin dynamics simulation using 
rather short chains 0|45| seem to support the linear be- 
havior. However, the above theories cannot explain the 
recent experimental result, namely that r ~ N 1,27 for 
polymer translocation through the solid-state nanopore 
|17|. Recently, Kantor and Kardar [3j| demonstrated 
that the assumption of equilibrium in Brownian polymer 
dynamics by Sung and Park 0] and Muthukumar |2^| 
breaks down more easily in the presence of a driving field 
and provided a lower bound N 1+ " for the translocation 
time by comparison to the unimpeded motion of the poly- 
mer. However, in their simulations they failed to observe 
this scaling behavior, which was attributed to the finite 
size effect of the polymer length. In addition, they 
also checked a polymer that is being pulled by one end. 
A quadratic dependence of t on N is predicted in their 
theory and confirmed by their simulations |33|. 

Most recently, we have investigated both free [3^| and 
forced (3(| translocation using the two-dimensional fluc- 
tuating bond (FB) model with single-segment Monte 
Carlo moves. For the free translocation, to overcome the 
entropic barrier without artificial restrictions we consid- 
ered a polymer, which is initially placed in the middle of 
the pore, and examined the escape time r e required for 
the polymer to completely exit the pore on either end. 
We found numerically that r e scales with the chain length 



N as r e ~ N 1+2v . This is the same scaling as predicted 
for the translocation time of a polymer which passes 
through the nanopore in one direction only EH • In 
addition, we also investigated the interplay between the 
pore length L and the radius of gyration R g . For L <C R g , 
we numerically verified that asymptotically r e ~ A 1+2ly , 
while for L 3> R g , we found r e ~ N. For forced translo- 
cation, we investigated the translocation dynamics under 
an external field within the pore. As our main result, we 
found a crossover scaling for the translocation time r 
with the chain length from r ~ N 2v for relatively short 
polymers to r ~ N 1+l/ for longer chains. Our results 
disagree with the experimental data that r depends lin- 
early on N in the case of a-hemolysin 0, EJ j but the 
predicted short chain exponent 1v sa 1.18 in three dimen- 
sions agrees reasonably well with the solid-state nanopore 
experiments of Storm et al. [nj , where the experimental 
setup is closer to the theoretical models investigated. No 
crossover in scaling behavior was observed in this experi- 
ment, presumably due to the fact that the polymers used 
in experiments are not long enough. 

In this paper we investigate the translocation dynam- 
ics using Langevin dynamics simulations (LP). Previous 
LD simulations using rather short chains |43l,l4fij seem to 
support the linear behavior for forced translocation time 
as a function of N, which is not expected to be correct for 
self-avoiding polymer. Also this result is in contradiction 
with the FB model simulations. (3(| On the other hand, 
the FB model is based on single-segment Monte Carlo 
moves and thus neglects all translational degrees of free- 
dom of the chain. Thus it is important to check that 
this does not adversely affect translocation dynamics. In 
section [n] we briefly describe our model and the simula- 
tion technique. In section ITTT1 we present our results for 
both free and driven translocation. For the free case, the 
emphasis is on verifying the scalin g ex ponent obtained 
previously with the FB-model [33L l35j. For transloca- 
tion under an electric field, we concentrate on studying 
the crossover in the scaling of translocation time as a 
function of the polymer length and on the influence of 
friction in the scaling exponent. Finally, the conclusions 
and discussion are in section HVI 



II. MODEL AND METHOD 

In the simulations, the polymer chains are modeled as 
bead-spring chains of Lennard- Jones (LJ) particles with 
the Finite Extension Nonlinear Elastic (FENE) potential. 
Both excluded volume and van der Waals interactions 
between beads are modeled by a repulsive LJ potential 
between all bead pairs: 



Ulj(t) 




< 2VV 
> 2 1/6 cr, 



(1) 



where a is the diameter of a bead, and e is the parameter 
adjusting the the depth of the potential. The connectiv- 
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ity between the beads is modeled as a FENE spring 

U FE N E {r) = ~\kR 2 ln(l - r 2 /R 2 ), (2) 

where r is the separation between consecutive beads, k 
is the spring constant and Rq is the maximum allowed 
separation between connected beads. 

In the Langevin dynamics method, each bead is sub- 
jected to conservative, frictional, and random forces Fj , 
Ff , and Ff , respectively, and obeys the following equa- 
tion of motion |4!| 

mr,=Ff + Ff + Ff , (3) 

where m is the monomer's mass. Excluded volume inter- 
action is explicitly included in Fp. Hydrodynamic drag 
is included through the frictional force, which for indi- 
vidual monomers is Ff = — £vi, where £ is the friction 
coefficient, and is the monomer's velocity. The Brow- 
nian motion of the monomer resulting from the random 
bombardment of solvent molecules is included through 
Ff and can be calculated using the fluctuation dissipa- 
tion theorem 

(Ff (*)> = 0; m 
(Ff (t) ■ Ff (f )> = 6k B T£S l3 5(t - t'). W 

The conservative force in the Langevin equation consists 
of several terms Ff = F LJ + F FE ne + F driving- Tnc 
driving force depends on potential difference. For the 
spontaneous case there is no driving force. 

The wall is described as I columns of stationary parti- 
cles within distance a from one another and they interact 
with the beads by the repulsive Lennard-Jones potential. 
Wall particle positions are not changed in the simula- 
tions. The pore is introduced in the wall by simply re- 
moving w beads from each column. Under an electric 
field, the potential is expressed as 

f -El, * > f/2; 

U e = I -Ex, 1/2 >x>-lh and y 2 < (w/2) 2 ; 
i E i x< -1/2, 

(5) 

where E is the electric field strength, I is the pore length 
and w is the pore width. 

In our simulations, the LJ parameters e and a fix the 
system energy and length units, respectively. Time scale 
is given by t^j = {mo 2 / e) 1 / 2 . The parameters are a = 1, 
Rq = 2a, k = 7e. The Langevin equation is integrated 
in time by a method described by Ermak and Buckholtz 
m| in 2D. For the pore, we set w = 2 and I = 1 un- 
less otherwise stated. To create the initial configuration, 
the first monomer of the chain is placed in the entrance 
of the pore. The polymer is then let to relax to obtain 
an equilibrium configuration such that the first monomer 
position is fixed at the entrance but the other monomers 
are under thermal collisions described by the Langevin 
thermostat. In all of our simulations we did a number of 
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FIG. 1: The distribution of 1000 escape times for a chain of 
length N/2 — 80 normalized by the mean value for free case. 
In the simulation, £ = 0.7 and fcsT = 1.2e. 

runs with uncorrelated initial states and random numbers 
describing the random collisions. The estimate for the 
translocation time was obtained by neglecting any failed 
translocation and then calculating the average duration 
of the accepted translocations. We note here that in some 
articles e.g. in Ref. |38| the estimate for the transloca- 
tion time is chosen to be the most probable translocation 
time. We checked that this does not change the scaling 
behavior of the chains. 



III. RESULTS AND DISCUSSION 

A. Free translocation 

The translocation of a polymer through a nanopore 
faces a large entropic barrier due to the loss of a great 
number of available configurations. In previous simu- 
lations [32l l33| . an artificial restriction is imposed to 
prevent the first monomer from escaping from the pore. 
However, this restriction is unphysical. Very recently, we 
suggested a new method to investigate the translocation 
dynamics without any restriction |35j. In this method, 
the middle monomer is initially placed in the middle of 
the pore. The polymer can escape the pore from cither 
side in time defined as the escape time r e . Here, we 
use this method to investigate the distribution of escape 
times and the scaling of r e with the chain length. As dis- 
cussed below, the distribution of escape times is similar 
to the distribution of translocation time in the study by 
Chuang et al. [32 • More importantly, escape time scales 
in the same way with the chain length as translocation 
time. 

In the absence of the driving force translocation is ex- 
tremely slow. For this reason we studied polymers of 
length only up to N = 300. The distribution of escape 
times is shown in Fig. ^for a polymer of length N/2 = 80 
normalized by its mean value. We find that the escape 
times are distributed on a wide range with a long tail. 
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FIG. 2: (a) Average escape time r e and (b) the most proba- 
ble escape time r as a fuction of the polymer length for free 
case with £ = 0.7 and ksT — 1.2e. The averages where cal- 
culated over 1000 averages for chain lengths up to N/2 — 80. 
For chains of length N/2 = 100 and 150 we had 300 and 100 
averages, respectively. 




FIG. 3: The polymer configurations for N = 300 under an 
electric field E = 5s, £ = 0.7 and ksT — 1.2e before translo- 
cation (a) and during translocation (b). 
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FIG. 4: The distribution of 1000 translocation times for a 
chain of length N = 100 normalized by the mean value in the 
case of an electric field of strength E — 5e as a driving force. 
The other parameters are £ = 0.7 and ksT — 1.2e. 

not all simulation runs lead to a successful translocation 
and even when they do, translocation times vary over a 
wide range of values. The success ratio describing the 
percentage of succesfull translocations depends on vari- 
ous factors. The driving force strength is the most im- 
portant factor to increase the success probability. On 
the contrary, increasing the chain length and the fric- 
tional force, decreases the success probability. Below, we 
consider only the translocation time r over all successful 
runs. 



The distribution is similar to those found previously with 
the FB model by Chuang et al. [J^] • The average escape 
time T e as a function of N is shown in Fig. . We find 
the scaling exponent x — 2.48 ± 0.07. It is clear that 
average escape time is different from the most probable 
escape time. However, we have checked that the scaling 
exponents are almost the same. Using the most probable 
escape time, we find x = 2.5 ± 0.1, see Fig. |2Jd). Above 
results are in excellent agreement with the theoretical 
prediction t ~ N 1+2v [32], and our previous results 

Te „ jy2.60i0.01 

B. Polymer translocation under an electric field 

For driven translocation, the chain is initially placed 
on one side of the pore with one end of it in the pore en- 
trance, as shown in Fig. Then the chain is allowed 
to reach an equilibrium state, but with the constraint 
that the first monomer is fixed. Once the polymer is in 
its equilibrium state, the first monomer at the entrance 
of the pore is released, and that moment is designated 
as the start of the translocation. The translocation time 
is defined as the time duration taken for the chain to 
move through the pore. The typical snapshots of poly- 
mer configurations during translocation under an electric 
field is presented in Fig. ^jp). It is important to note that 



1. Distribution of translocation times 

As shown in Fig. 0] the distribution of translocation 
times has a qualitatively different shape compared with 
the distribution of the escape times in Fig. ^ an d the 
the distribution of translocation times studied by chuang 
et al. 32] for the free case. In this case the distribution 
of translocation times is symmetric and narrow without 
a long tail. The stronger the diving force, the narrower 
distribution of translocation times. A very important 
feature is that the average translocation time and the 
most probable translocation time are almost identical. 
Thus the average translocation time r is well defined. 



2. Effects of £ and E on t 

In the Langevin equation, the friction coefficient £ de- 
scribes the strength of coupling to the solvent. The value 
£ = corresponds to the absence of the solvent molecules 
and would result in ballistic movement of particles. In 
the opposite limit £ — > oc, i.e. the overdamped limit, 
inertia plays no role in the dynamics. 

The translocation time as a function of friction £ for 
E = 5e and ksT — 1.2e is shown in Fig. [5J As expected, 
r increases with increasing £ due to the increases of the 
frictional force. An important feature is that two regimes 
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FIG. 5: The translocation time as a function of the friction 
for a chain of length N = 100 under an electric field E = 5e 
and k B T = 1.2e. r is an average of 1000 translocations. The 
slope of the straight line in the figure is 1. 



are observed. In the first rigime, when £ < 0.5, the 
translocation time increases very slowly with £. However, 
in the second regime, for £ > 0.5, it increases rapidly and 
linearly with £. These results can be understood from 
the interplay between the driving force and the frictional 
force according to the Langevin equation. For small £, 
the driving force and inertia play dominant roles com- 
pared with the frictional force. As a result, the translo- 
cation time t has a weak dependence on friction, thus 
increases only slightly with £. However, when £ is large 
enough, inertia can be neglected. According to the bal- 
ance of the driving force and the frictional force, the hor- 
izontal velocity of the center of mass v cmx ~ y. Thus 

the translocation time r ~ — — 



I l5£ 



The influence of the driving force on the translocation 
time has also been investigated and the results are illus- 
trated in Fig. El for a chain of length TV = 100. We 
find that r ~ £?-°- 97 ± - 02 which is in excellent agreement 
with simple argument above that r ~ E~ x in the regime 
where inertia can be neglected. This finding is in excel- 
lent agreement with the results from the experiments of 
Kasianowicz et al. 0, who found that the translocation 
is inversely proportional to the electric field strength. By 
contrast, our result disagrees with another experimental 
finding of an invers e q uadratic dependence of E for a- 
hemolysin channel 



so f|l 

m 



5 
4.5 
4 

•3.5 
' 3 

2.5 
2 

1.5 



4000 r 

} 

f 2000 



10 

E / e 



20 



x=-0.97±0.02 



0.5 1 
log10(E)/ E 



1.5 



FIG. 6: The translocation time as a function of electric field 
strength for a chain of length N — 100. r is an average of 
1000 averages. Here k B T = 1.2e and £ = 0.7. 
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FIG. 7: The translocation time as a function of polymer 
length N for £ = 0.7, (a) E = 5e, and k B T = 1.2e, (b) 
E = 2.4e, and k B T = 1.2e, (c) E = 8.3e, and k B T = 2e, (d) 
E = 5e, k B T = 1.2e, and w = 3. r is an average of 1000 runs. 
In each figure, a shifted solid line is plotted beneath the data 
to which the curve is fitted. The curve is continued with the 
same slope but for clarity it is plotted as a dashed line. The 
inset shows the local scaling exponent x in which each point 
is the slope of a curve fitted to four points. 



3. Crossover scaling behavior of the translocation time with 
the chain length 

The translocation time as a function of the polymer 
length with £ = 0.7 is presented in Fig. [7\ In each figure, 
a shifted solid line is plotted beneath the data to which 
the curve is fitted. The curve is continued with the same 
slope but for clarity it is plotted as a dashed line. The 
inset shows the local scaling exponent x in which each 
point is the slope of a curve fitted to four points. 

The result for E — he and /c^T = 1.2e is shown in Fig. 



EH)- One of the main features is that a crossover scaling 
behavior is observed. For short chains (N < 70), the scal- 
ing exponent x = 1.50 ±0.01, which is in good agreement 
with 2v, where v = 0.75 is the Flory exponent for a self- 
avoiding walk in 2D. However, for N > 300 the exponent 
becomes x — 1.69 ± 0.04, which is close to 1 + v = 1.75. 
The polymer assumes a more compact configuration re- 
sulting from both the dynamical effect and geometrical 
restriction. This leads to a smaller effective value for the 
exponent v as compared to its Flory value • For the 
different electric strength E and temperature T and pore 



6 



x 1.6^***" 

1.4-1 ■ — 

500 
N 



(a) 



log10(N) 

(b) 



(o) 




_ - x=1 

1.8 


67±0.03 


x 1.6 






x=1.64 ±0.01 , , 






' 500 




N 



loglO(N) (C) 

FIG. 8: The translocation time as a function of polymer 
length N with an electric field of strength E — 5e as a driving 
force for k B T = 1.2e and (a) £ = 0.4, (b) £ = 1.5 and (c) 
£ = 3. r is an average of 1000 runs. 



width w, the same exponent and crossover behavior is 
observed, as shown in Fig. d>), c) and d), respectively. 
Next, we investigate the effect of the friction on the scal- 
ing exponent. In Fig. 00 we present results for r as a 
function of N with E — 5e for three different frictions: 
£ = 0.4, £ = 1.5 and £ = 3. For £ = 0.4 and £ = 1.5, a 
clear crossover is observed, but for £ = 3 only r ~ 7V 1+2!y 
occurs. 

To understand the above results, one should consider 
the configuration of the monomers after exiting from the 
during pore under an electric field, as shown in Fig. Et>). 
For driven translocation, there is a qualitative difference 
compared with the purely diffusive free translocation for 
long chains: translocation under a strong driving force 
is much faster than that for free diffusion and thus the 
translocated monomers don't have time to diffuse away 
from the vicinity of the pore exit, see Fig. |3Jd). As 
to the effect of the friction on the translocation dynam- 
ics, on the one hand, increasing £ leads to the increase 
of the frictional force, resulting in slowing down of the 
translocation. On the other hand, the Rouse relaxation 
time for the self-avoiding chain is proportional to jr^fi 
which rapidly increases with £. The latter factor is al- 
ways dominant and thus the translocation time is much 
shorter than the Rouse relaxation time. As a result, it is 
more difficult for translocated monomers to diffuse away 
from the pore exit even for short chains with higher £. 
We checked that the radius of gyration is larger before 
translocation than that immediately after it. This fact in- 
dicates the higher density of translocated monomers near 
the pore exit. It can be imagined that the higher den- 
sity of translocated monomers near the pore exit greatly 
slows down the translocation. Therefore, for intermedi- 
ate values of £, two scaling regimes for short and long 



chains are observed, respectively. However, for high fric- 
tion such as £ = 3, the Rouse relaxation time is so long 
compared to the translocation time that it is more diffi- 
cult for the translocated monomers to diffuse away even 
for short chains, thus the translocation dynamics directly 
enters into the regime where r ~ N 1+2v . 

Our results disagree with the experimental data that 
r depends linearly on N for translocation through a- 
hemolysin pore Q, , but the predicted short chain ex- 
ponent 2v = 1.18 in 3D agrees well with the solid-state 
nanopore experiments of Storm et al. 17], who found an 
exponent of 1.27. As pointed out by Storm et al. [l7| . 
the linear behavior for cv-hemolysin pore may be from 
the specific and complicated interaction between DNA 
chain and the pore. The beginning of the crossover re- 
gion occurs at N > 300 which is beyond or near the 
longest ssDNA and dsDNA used in the experiments so 
far Thus, it is not surprising that crossover 

in scaling behavior has not been experimentally observed 
yet. Theoretically, as expected our results disagree with 
the linear dependence r ~ N prediction by Sung and 
Park |19| and Muthukumar |22j| . which is invalid for 
nonequilibrium translocation. However, our results also 
disagree with the previous Langevin dynamics simula- 
tions for relatively short polymers pM 0] , which show 
r ~ N . The present results agree with the main find- 
ings in our previous FB-studies 36] for both short and 
long chains and the theoretical prediction r ~ N 1+l/ by 
Kantor and Kardar 1331. 



4- Crossover behavior of the translocation velocity 

To study translocation dynamics in detail, we also cal- 
culated the translocation velocity v as a function of poly- 
mer length N. The translocation velocity can be mea- 
sured in several ways. A simple way is to measure the 
average horizontal velocity of the center of mass of the 
polymer over the whole duration of all successful runs 




(6) 



where v xi is the horizontal component of the velocity of 
monomer i. In addition, the translocation velocity can 
also be defined as 



(go) 

in)' 



(7) 



where Xq is the horizontal coordinate of the last monomer 
in the initial equilibrated configuration of the polymer 
before the translocation and n is the translocation time 
for every successful run. 

According to our numerical results for crossover scaling 
behavior, we must have that 
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FIG. 9: The translocation velocity as a function of polymer 
length N with an electric field of strength E — 5e, ksT — 1.2e 
and (a) £ = 0.7 and (b) £ = 3. Averages are taken over 100 
runs, xl corresponds to the velocities defined as in Eq.0 
and the x2 to velcities defined as in Eq.JHJ. The latter results 
have been shifted for clarity. 



N = 300 are presented in Fig. ^| 

The waiting time depends strongly on the monomer 
position in the chain. For the shorter polymer N = 100 
the monomers approximately in the middle of the poly- 
mer need the longest time to translocate and the distri- 
bution is close to symmetric. However, for the polymer of 
length N = 300, it's approximately the 220^ 1 monomer 
that needs the longest time to thread the pore, which in- 
dicates that during late stages of translocation the high 
density of segments of a long polymer near the pore 
slows down the translocation. The same phenomenon 
was found in our earlier study with the FB-model |3r1 |. 



IV. CONCLUSIONS 



monomer number 



(a) 



(b) 



FIG. 10: The waiting times for a polymer of (a) length N = 
100, (b) N = 300 for E = 5s, k B T = 1.2e and £ = 0.7. The 
waiting times are averages of 1000 runs. 



In Fig. Et), we present the polymer velocity as a func- 
tion of chain length for E = 5e, ksT — 1.2e and £ = 0.7. 
We get v ~ jV-°- 778±a003 and v ~ A- a842±a03 for short 
chains and v ~ Ar-i-02±o.05 and v ~ AT-i-00±o.04 for long 

chains, respectively. There is a clear crossover in the 
translocation velocity from v ~ N~ v to v ~ N , which 
is in agreement with Eq.(|HJ. This simple test confirms 
that the crossover in the translocation time takes place 
because of a crossover in the translocation velocity. Here, 
we should mention that for high friction such as £ = 3, 
there is no crossover and only v ~ N^ 1 is observed, see 

Fig. BO . 



C. Waiting times 

The dynamics of a single-segment passing through the 
pore during translocation is an important issue, as far as 
experiments are concerned. The non-equilibrium nature 
of translocation has a considerable effect on this. We 
have numerically calculated the waiting times for each 
monomer passing through the pore. We define the wait- 
ing time of monomer s as the average time between the 
events that monomer s and monomer s + 1 exit the pore. 
The waiting times for chains of length N = 100 and 



In this work, we investigate the dynamics of poly- 
mer translocation through a nanopore using the two- 
dimensional (2D) Langevin dynamics simulations. For 
free translocation, we consider a polymer which is ini- 
tially placed in the middle of the pore and study the 
escape time r e . The distribution of the escape times is 
wide and has a long tail. We verified that average es- 
cape time (as well as the most probable escape time) 
scales with the chain length N as r e ~ N 1+2u . These re- 
sults confirm previous theoretical studies based on a less 
miscrosopic fluctuating bond model. For forced translo- 
cation, we concentrate on studying the influence of the 
friction £, driving force E and polymer chain length N . 
For strong driving forces, the distribution of transloca- 
tion times is narrow without a long tail and symmetric, 
which is completely different compared with free translo- 
cation. The influence of the £ on the translocation time 
depends on the ratio between driving force and frictional 
force. We find translocation time r is inverse propor- 
tional to driving force for a wide range of friction values. 
Finally, as regards to the dependence of translocation on 
the length of the polymer, we find a crossover scaling for 
the r with N from r ~ N 2v for relatively short poly- 
mers to r ~ N 1+l/ for longer chains at moderate values 
of the friction. For higher £, there is no crossover be- 
cause in this case the Rouse relaxation time is extremely 
long compared with translocation time. These scaling 
behaviors can be understood from the observation that 
in the limit of large N and £, there is a high density of 
segments near the exit of the pore, which slows down the 
translocation process due to slow relaxation of the chain. 
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